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Abstract. For certain hierarchical structures, one can study the percolation problem 
using the renormalization-group method in a very precise way. We show that the idea 
can be also applied to two-dimensional planar lattices by regarding them as hierarchical 
structures. Either a lower bound or an exact critical probability can be obtained with 
this method and the correlation-length critical exponent is approximately estimated 
as v s» 1. 
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The percolation problem is a question about how a global connection can be made 
possible by randomly filling local components by a certain probability p. While it 
can be explained in purely geometric terms without any interaction, when a global 
connection actually appears, the macroscopic behavior of the system exhibits all the 
characteristic features of a continuous phase transition with a diverging correlation 
length, just as we observe in other interacting spin systems such as the two-dimensional 
(2D) Ising model pQ. This analogy is given a precise meaning by the Fortuin-Kasteleyn 
representation of the g-state Potts model |2j, where the percolation turns out to be 
equivalent to the limit of q — > 1. Since the percolation transition at a critical probability 
p c has a diverging correlation length, every microscopic length scale becomes irrelevant 
with respect to the critical phenomena, and the system behaves as if it does not have any 
specific length scale. This is a qualitative explanation of the reason why a percolating 
cluster connecting two opposite sides of a 2D plane has a fractal dimension at p = p c . 
The lack of a specific length scale implies that the system remains statistically invariant 
even if we zoom the system up or down, and this scale invariance readily lends itself to 
a renormalization-group (RG) study of the percolation problem [31 HJ [5]. 

In certain cases where the underlying structure itself is fractal, it is possible to 
carry out the RG calculation to a good approximation or exactly, exploiting this fractal 
property [6j [7]. Such fractal structures usually contain groups of bonds which connect 
longer and longer distances in a regular fashion. For this reason, one can sometimes 
arrange the groups of bonds in a hierarchical way according to their connection lengths. 
Figure [H(a) is an example of a hierarchical structure called the enhanced binary tree, 
which is obtained by adding horizontal bonds to the simple binary tree. It is hierarchical 
in the sense that filling a horizontal bond is comparable to a very long connection along 
the bottom layer and the connection length is dependent on the level of the horizontal 
bond [8]. That is, a horizontal bond in the highest level can connect two points at 
distance 7 along the bottom layer at maximum. For a horizontal bond at the next 
highest level, this maximum connection distance is only as large as 3 lattice spacings. 
An RG scheme for the enhanced binary tree is described in [8] as shown in figure E^b): we 
calculate the probability for any of the leftmost points to connect to any of the rightmost 
points within the cell as a function of the bare coupling p and a coarse-grained effective 
coupling z n , and then replace this probability by a new effective bond with strength 
z n+ i. The resulting expression for z n+1 is written as 

z n+1 = p + (1 - p) [(1 - pfzl + 2p(l - p)z 2 n + p 2 z n ] . 

By asking when z n = z n+ \ = z^ becomes 1, we obtained a lower bound of the percolation 
threshold as p c > 1/2 [8J, which is consistent with the conclusion in j9] that p c = 1/2. 
Note that we get a lower bound since in iterating z n to z n+ i, there is a small chance to 
regard a layer as percolated when it is actually not [see, e.g., figure Hfc)], whereas the 
opposite is not possible. 

Although the above RG scheme is devised to investigate a hierarchical structure, 
we show in this work that it can be applied to non-hierarchical planar lattices as well. 
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Figure 1. (a) Enhanced binary tree, a hierarchical structure derived from a simple 
tree with branching number 2. (b) RG scheme for the enhanced binary tree, where the 
connectivity over the cell above is coarse-grained into a single bond hlled by probability 
z n+ i. (c) This layer is not connected from left to right even though the cells inside it 
appear as filled according to the recursion scheme in (b). The solid and dotted lines 
represent filled and empty bonds, respectively. 



In figure [2^a), we present a variation of the RG scheme shown above. The similarity 
is obvious: we have taken away only one bond out of those in figure [Ob), and this is 
meant to describe the triangular lattice. It leads us to the following recursion, 

Z n +1 =P+ (1 ~p) \p+ (1 -P)z n ] 2 . (1) 

Again, the bond connection in lower levels, composed of p and z n , is converted to a 
single bond with z n+ i at a higher level. This distinction of levels might look arbitrary 
since the bonds in the plane do not have any hierarchy. However, the important point 
is that all the argument above to find a lower bound remains still legitimate from this 
viewpoint. Solving ([[]) for z n+ i = z n = z^, we find that 

_ p(l +p-p 2 ) 

Zco ~ (i- P r 

and consequently, — 1 at p* — 1 — l/y/2 w 0.293. Comparing this to the exact 
bond-percolation threshold in the triangular lattice, p\ ~ 0.347 [TO], we see that our 
method indeed yields a lower bound. We now extend the cells to be renormalized by 
adding one more level. That is, let us denote the width of the cell as w and consider 
the case of w — 2. For the triangular lattice, the shape of such a larger cell is given in 
figure EJjb). By enumerating all the possible cases, the recursion relation is obtained as 

z n+1 = 3p 9 zl - 25/4 + 90p 7 4 - 182p 6 4 

+ 224p 5 4 - 168p 4 4 + 70p 3 4 - I0p 2 z 3 n 

- 3pz 3 n + zl- 7p 9 z 2 n + 53p 8 z 2 n 

- 171p 7 z 2 n + 303p 6 z 2 n - 315p 5 z 2 n + I87p 4 z 2 n 
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Figure 2. (a) A variation of the RG scheme in figure QJb). Here it describes the 
triangular lattice with width w — 1. (b) A larger cell with w — 2 for the triangular 
lattice. The RG scheme for the honeycomb lattice with (c) w = 1 and (d) w = 2 can be 
constructed in the same way, as well as that for the square lattice with (e) w = 1 and 
(f) w = 2. The double lines represent coarse-grained effective bonds and the dotted 
lines in (c) to (f ) mean that the connections do not correspond to any bare interaction. 



- 53p 3 4 + P 2 z 2 n + 2pz 2 n + 5p 9 z n 

- 33p s z n + 89p 7 z n - 121p 6 z n + 79p 5 z n 

- llp A z n - I3p 3 z n + 5p 2 z n - p 9 + 5p s 

- 8p 7 + 12p 5 - 8p 4 - Ap 3 + 4p 2 + p, 

and we find its limiting value as 

with Fi{p) = 4 P 6 - !6p 5 + 21p 4 - 6p 3 - 6p 2 + 2p + 1, F 2 (p) = Ap 12 - 40p n + 
176p 10 - 400p 9 + 653p 8 - 508p 7 + 48p 6 + 236p 5 - 126p 4 - 32p 3 + 20p 2 + 8p + 1, and 
F 3 (p) = 6p 6 -32p 5 +66p 4 -64p 3 + 26p 2 -2. The solution of = 1 is found at p* w 0.300, 
which is an improved lower bound compared to the previous one, p = 1 - 1/V2 w 0.293, 
even though the convergence turns out to be rather slow. 

If we also regard the honeycomb lattice as hierarchical, we can consider an RG 
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scheme as depicted in figure El^c). By calculating the probabilty for any of the leftmost 
points to connect to any of the rightmost points within the cell, we find 

z n+l = [p+ (I -p)Z n f . (2) 

By a little algebra as above, we find z^ = p 2 /{l — p) 2 , which becomes one at p* = 1/2. 
Again, this is lower than the exact value p\ ~ 0.653 [10]. One may expect an improved 
estimate by considering a larger cell shown in figure EJ^d), which leads to 

z n+l = - Ap 7 z 3 n + 17/4 - 26p 5 4 + I5p 4 z 3 n 

- p 2 z 3 n - 2pz 3 n + z\ 

+ 10p 7 zl - 35p 6 zl + A0p 5 z 2 n - Y2p 4 z 2 n 

- Ap 3 z 2 n - p 2 z 2 + 2pz 2 - 8p 7 z n 

+ 21p 6 z n - 12p 5 z n - Qp 4 z n + Ap 3 z n 
+ p 2 z n + 2p 7 - 3p 6 - 2p 5 + 3p 4 + p 2 . 

The limiting solution is 

gi(p) ~ VG^pj 
G 3 ( P ) 

where G x {p) = -6p 5 + Qp 4 + Ap 3 - p 2 - 2p - 1, G 2 (p) = V° - lQp 9 + 2Ap 8 - 12p 7 + 
12p 6 - 20p 5 - lip 4 + Ap 3 + 10p 2 + Ap + 1, and G 3 (p) = 8p 5 - 18p 4 + 8p 3 + Ap 2 -2. We find 
that Zoo = 1 at p* ~ 0.537. Using the duality relation p\ +p^ = 1 [ID] , we may turn this 
result to an upper bound of the bond-percolation threshold for the triangular lattice. 
That is, our method gives a possible region of the threshold as 0.300 < p\. < 0.463, or 
equivalents, 0.573 < p\ < 0.700. 

A more interesting case is found by considering the horizontal bonds in figure E^a) 
and figure [2](b) as fictitious [figure [2]^e) and figure Eff)]- This corresponds to the square 
lattice, and the interaction in the horizontal direction will appear only as an effective 
one mediated by shorter bonds. Then we can simplify (OQ) as 

z n+1 = [p+ (1 -p)z n } 2 , 

which happens to be the same as ([2]). Therefore, we find p* = 1/2 once again, but 
this value is identical to the exact value for the bond-percolation problem in the square 
lattice [11]. Since this method is supposed to give a lower bound, it should not be 
possible to improve this result further, so it will be worth checking whether this value 
really remains unchanged for a larger cell. From a larger cell depicted in figure [2](f), we 
obtain a recursion 

z n+l = - 3p 6 z 3 n + Up 5 z 3 - 25p 4 z 3 + 20p 3 z 3 

- 5p 2 z 3 n - 2pz 3 n + z 3 n + 7p 6 z 2 n 

- 28p 5 z 2 + A0p 4 z 2 n - 22p 3 z 2 n + p 2 z 2 n 
+ 2pz 2 n - hp Q z n + 16p 5 z n - \Ap 4 z n 
+ 3p 2 z n + p 6 - 2p 5 - p 4 + 2p 3 + p 2 , 
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and find its limiting value as 

Hi(p) - y/Hjp) 
H 3 (p) 

with H x {p) = 4p 4 -6p 3 -p 2 +2p+l, H 2 (p) = Ap 8 -lQp 7 + 27p 6 -12p 5 -I5p 4 + Qp 2 +Ap+1, 
and H 3 (p) = 6p A — 16p 3 + 12p 2 — 2. The critical value making — 1 is also = 1/2, 
as expected. The fact that p* does not change with w could be an evidence that the 
bond-percolation threshold is located exactly at p = 1/2 for the square lattice. 

In addition, we can argue that the connection probability over distance I would 
be roughly determined by (z^) = e ll ° sz °° near the critical point. In other words, 
the correlation length would be written as £ ~ — 1/ logz^. The slope of z^ around 
p = p c does not vanish in every case considered above, so it generally behaves as 
Zoo ~ a(p — Pc) + 1 where a = dz^/dpl » ~ 0(1) at p = p c — e with positive e C 1. 
Therefore, we see that 



log Zoo log [a (p -p c ) + 1] 

~ (Pc - P)"\ 

by using log(l — ae) ~ — ae. Since the correlation length is assumed to diverge as 
£ ~ |p— p c \~ u , this argument gives us an approximate estimate of the critical exponent as 
v ~ 1, which is an underestimate compared to the exact value, v = 4/3 [12J. It is worth 
noting that this RG scheme does not make use of any explicit scaling transformation: 
we do not zoom up or zoom down the system at criticality as usually found in RG 
studies [31 HI E]. In arguing the value of u, therefore, we evaluate it directly in units of 
the given lattice spacing instead of any zooming ratio. By setting z n = z n+ i, in a sense, 
it is the translational invariance that we are actually exploiting in this study. 

In summary, we have shown that the RG scheme devised for a hierarchical structure 
can be also applied to the 2D lattices even though they are not hierarchical. It generally 
yields a lower bound, but correctly predicts the bond-percolation threshold for square 
lattice. We have also approximately estimated v ~ 1. This method is more related to 
the translational invariance rather than to the scaling invariance at criticality. 
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